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Abstract 

Various phase concepts may be treated as special cases of 
the maximum Hkehhood estimation. For example the dis- 
crete Fourier estimation that actually coincides with the op- 
erational phase of Noh, Fougeres and Mandel is obtained for 
continuous Gaussian signals with phase modulated mean. 
Since signals in quantum theory are discrete, a predic- 
tion different from that given by the Gaussian hypothesis 
should be obtained as the best fit assuming a discrete Pois- 
sonian statistics of the signal. Although the Gaussian es- 
timation gives a satisfactory approximation for fitting the 
phase distribution of almost any state the optimal phase 
estimation offers in certain cases a measurable better per- 
formance. This has been demonstrated in neutron-optical 
experiment. 



1 Introduction 

Physics enables us to comprehend Nature by considering 
intimate relations between various effects. Any physical 
observation can always be compared and analyzed in re- 
lation with a particular internal model, providing us with 
some additional insight into the laws of Nature. However, 
this effort need not and usually does not tend to a unique 
solution. It may happen that there are several plausible 
models and the given observation is not able to discrimi- 
nate among them. On the other hand it may also happen 
that some of the assumptions about the system need not 
to be apparent. Certain statements therefore pretend to be 
more general than they are in reality. This interplay be- 
tween physics and philosophy may be demonstrated on the 
long standing problem of quantum theory-on the problem 
of quantum phase. Phase measurements belong to standard 
detection techniques revealing the wave property of the sig- 
nal. The quantum phase, however, encountered theoretical 
difhculties when an adequate quantum theory has been con- 
structed 1^. 



There are several concepts for the description of phase 
in quantum theory at present. For an up to date overview 
see § |. Some of them are accenting the theoretical as- 
pects, other the experimental ones. The operational ap- 
proach formulated by Noh, Fougeres and Mandel (NFM) 
[|[ 1^ is motivated by the correspondence principle in classi- 
cal wave theory. The interference pattern is adopted for the 
scheme where the sin and cos function of the phase shift are 
measured simultaneously in the 8-port homodyne detection. 
An equivalent measurement may be realized on the Mach- 
Zehnder interferometer, provided that the measurement of 
an unknown phase shift is done with and without an ad- 
ditional 7r/2 phase shifter. The NFM scheme is plausible 
whenever the signal behaves like a classical wave since be- 
sides the principle of correspondence, no other assumptions 
has been used. 

As a particular result, the NFM scheme provides the op- 
timum result, provided that the statistics of the signals is 
represented by Gaussian statistics with a phase sensitive 
mean and a phase insensitive noise. Since realistic signals 
in the quantum world do not meet these rather restricting 
criteria, phase prediction based on them is not optimum in 
general. The difference between Gaussian estimation and 
optimum treatments is caused by the statistical nature of 
the phenomena. This can be experimentally registered in 
an interferometer with discrete Poissonian signal. Several 
possible interpretations of this comparison are noteworthy: 
i) This may be considered as a testing of the operational 
quantum phase prediction. It quantifies how well the NFM 
phase concept fits reality, ii) It may be interpreted as a non- 
trivial statistically motivated "quantum calibration" of an 
interferometer. The visibility of interference fringes is usu- 
ally used for this purpose. However, this criterion focusses 
on the wave property of the detected signal only. The pro- 
posed method involves and evaluates the whole detection 
process, particularly the ability to control phase shift in 
the experimental arrangement and the statistics of the de- 
tected signal. This seems to be in accordance with the prag- 
matic interpretation of quantum theory, where the results 
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depend irreducibly on both state preparation and measure- 
ment, iii) In the framework of wave-particle duahty, the 
proposed treatment tries to answer the question: "Does the 
interfering signal resemble more discrete particles or classi- 
cal continuous waves ?" iv) It provides an example of in- 
direct observation of several parameters. Particularly, by 
detecting the interference fringes, the phase shift as well as 
the visibility may be determined as fluctuating variables. 
It provides one of the simplest examples of the so called 
"quantum state reconstruction" procedure. 

This paper is organized as follows. In the second sec- 
tion, the main idea of comparing the NFM scheme with an 
optimum phase prediction is developed. In the third sec- 
tion, the idea is generalized in order to apply the scheme 
for phase measurement in matter wave interferometry. The 
fourth section provides the experimental realization in neu- 
tron interferometry the results obtained. 

2 Statistical formulation of opera- 
tional phase concepts 

In this section the operational phase concept will be natu- 
rally embedded into the general scheme of quantum estima- 
tion theory |^ . A similar approach has been already used 
in ^, H]. However, since the purpose of the detection scheme 
is to predict the phase shift after each run, the point esti- 
mators of phase corresponding to the maximum-likelihood 
(MaxLik) estimation will be used here Assume an 

ideal device with four output channels enumerated by in- 
dices 3,4,5,6, where the actual values of intensities are reg- 
istered in each run. The values fluctuate in accordance with 
the statistics of continuous Gaussian signals. The mean in- 
tensities are modulated by a phase parameter 9 as 



h.i^-{l±Vcos9) 



(1) 



where / and V are total input intensity and visibility of 
the interference fringes, respectively. The energy is split- 
ted symmetrically between all the output ports. This de- 
vice represents nothing else than a classical wave picture 
of the original 8-port homodyne detection scheme. Equiv- 
alently, it also corresponds to a Mach-Zehnder interferom- 
eter, when the measurement is done for an unknown phase 
shift together with a zero and a 7r/2 auxiliary phase shifters. 
In this case, the data are not obtained simultaneously, but 
should be collected during repeated experiments. Provided 
that a particular combination of outputs /a , J4 , /s , /g has 
been registered, the phase shift should be inferred. In ac- 
cordance with the MaxLik approach ||l^, the sought-after 
phase shift is given by the value that maximizes the like- 
lihood function. The likelihood function corresponding to 
the detection of given data reads 
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[h-hf-ih-Ie]')}. (2) 



Here the variation a represents the phase insensitive noise 
of each channel. The sampling of intensities may serve for 
an estimation of phase shift, the average number of particles 
and the visibility simultaneously. A notation analogous to 
the definition of phase by Noh, Fougeres and Mandcl can 
be introduced as 



The likelihood function may be rewritten to the form 
C{9,VJ)^ exp{--ij[/-i^/,]2} 
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and is maximized by the choice of parameters 
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Hence the operational phase concept of Noh, Fougeres and 
Mandel is nothing but the MaxLik estimation for waves 
represented by continuous Gaussian signal with phase- 
independent and symmetrical noises. These rather strong 
assumptions are associated with the behaviour of waves in 
classical field theory. 

Since realistic signals are discrete they meet neither of 
these criteria and therefore, deviations in the optimum 
phase prediction should be expected. Assume the Poisso- 
nian statistics of an ideal laser. Together with the phase, all 
the parameters which are not controlled in the experiment 
will be optimally predicted as well. Denote for concreteness 
the detected discrete values as numbers 71.3,714,715,716. The 
likelihood function corresponding to this particular detec- 
tion as a function of the parameters 9, V and N reads 

ri3+n4+n5+»i6 
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The MaxLik estimation for parameters gives the optimum 
values for the phase shift, the visibility and the mean par- 
ticle number as 
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These relations provide a correction of the Gaussian wave 
theory with respect to the discrete signals. Besides the 
phase shift, the visibility of the interference fringes and the 
total energy input can be evaluated simultaneously. 

An apparent difference between relations (^-^) and 
|l2|) represents the theoretical background of the presented 
treatment. Adopting the interpretation of Ref. in both 
these approaches the unnormalized cos and sin functions of 
the phase shift are detected. However, the normalizations 
differ in both approaches. In the former Gaussian case, the 
normalization is performed only once, whereas in the latter 
Poissonian case it is done in two steps. The cos and sin func- 
tions of phase are normalized separately with respect to the 
total number of particles on both the output ports and then 
again among themselves. Obviously, both predictions will 
coincide provided that there is almost no information avail- 
able in the low field limit I,N ^ 0. Similarly in the strong 
field limit I, N cxd, where any statistics approaches the 
Gaussian one, the differences must disappear. Possible devi- 
ations may appear in the intermediate regime, characterized 
approximately by conditions I,N ^ 0(1). The test of the 
difference between (^) and ( |lO| ) is proposed as controlled 
phase measurement. The phase difference may be adjusted 
to a certain value and estimated independently using both 
the methods and ([l^) in repeated experiments. The ef- 
ficiency of both methods is then compared by evaluation of 
confidence intervals. Since any imperfections of the detec- 
tion scheme will smoothen the differences, it is questionable 
whether both schemes can be experimentally distinguished. 
This idea will be pursued in the following sections. 

Before doing this, the statistical analysis may clarify some 
subtle points of the NFM treatment, particularly then the 
nature of discarded data. Obviously the data yielding the 
ambiguous phase e*^"^" = 0/0 in equation (|^) would pro- 
vide zero visibility for both (j^) and (|ll]) . For a more detailed 
analysis, the Bayes theorem may be applied as well. In this 
case the likelihood functions quantifies the phase informa- 
tion involved in the detected data as posterior phase dis- 
tribution. Gaussian statistics provides homogeneous poste- 
rior phase distribution, whereas Poissonian statistics yields 
the four peak posterior distribution of the phase shift re- 
sembling effectively the homogeneous one. This statisti- 
cal analysis supports the conclusion of Refs. |]l3|, |l^, |l^ 
that the detected data cannot be discarded. Provided that 
some data are discarded, the average number of particles, 
i.e. the average energy corresponding to the phase detec- 
tion changes. Particularly, provided that the experiment 
has been done AI times and a total number Mi of parti- 
cles has been detected in each run, the average number of 
particles is simply J2i all -^i/^- Provided that some data 
appearing Aid times are discarded and not included in the 
evaluation the corresponding phase estimation is done with 
an average intensity loss Ci.aH ^i-Y.iMsc Mi) / {A/I - Aid) ■ 
Obviously, the measurement with and without discarding is 
done with a different energy input and this difference may 
be substantial. For example, this explains the ambiguity of 
the interpretation in Ref. [|6| . When the phase of a quan- 
tum field is measured against a classical field, no data are 



discarded since the field is strong. The relative phase of two 
quantum fields may then be evaluated separately against a 
strong field and a difference of two such phases values pro- 
vides the relative phase. However, this cannot be directly 
compared with direct measurements of two (weak) quantum 
fields against themselves, when some data are canceled. Al- 
though both measurements have been done with the same 
quantum states, the average energies of the corresponding 
phase detection may differ significantly. No wonder that 
the measurement with higher average energy gives better 
(sharper) results as observed in Ref. JlGf . 

The analysis given here will be adopted to the case of 
neutron interferometry, when the number of particles can be 
detected and discriminated with high efficiency (l^ . Highly 
efficient neutron detectors then provide almost perfect neu- 
tron number measurements. However the interference pat- 
tern exhibits lower visibility and the splitting process is far 
away from the symmetrical case. This will be taken into 
account when adopting the scheme to the case of neutron 
interferometry. 

Analogous analysis for detection of light needs other mod- 
ification. Since the efficiency of photodetection is less than 
unity, number of particles cannot be detected in such a 
straightforward way. Instead, the absence or presence of 
the week signal on the detector can be used for phase esti- 
mation. However, this will be dealt with separately in the 
forthcoming publication. 

3 Phase estimation in neutron in- 
terferometry 

Assume the following modification of the scheme developed 
above. The output ports of an interferometer are considered 
to be nonsymmetrical. The measurement is done with many 
auxiliary phase shifters. Neutron beams inside the perfect 
crystal neutron interferometer will be described at first as 
a classical wave. The signals 1°, /j* detected at the two 
output ports o,h (Fig. 1) are regarded as stochastic Gaussian 
intensities with phase sensitive means 

/; = 1° + cos{e + Aj) 

I'/ = I'' -I^ cos{9 + Aj), (13) 

where 6 is the true phase shift between the two branches of 
the interferometer, 1° and are the mean intensities of the 
two interference fringes, stands for the modulation am- 
plitude of the interference fringes (unnormalized visibility) 
and Aj are the values of the auxiliary shifts 

A, = j=O..N-l, NeM. (14) 

The Gaussian statistics of the detected signal yields its like- 
lihood function in the form 

ae) cx cxpj-— ^ [(/; - /;)^ + (i^ - i^)']j. (is) 
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Introducing the complex parameter R as 

R = ^(/;-/j')exp(-zA,), (16) 
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the parameters maximizing the Hkchhood function read 
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\ N ' ' 



(17) 
(18) 

(19) 



The relation (19) follows from the condition of semidefinite- 
ness of the amplitude < min{/°,/''}. Expression ( |l7| ) 
represents a generalization of the NFM formula (H), which 
may be recovered provided that the measurement is done 
at the two positions Aq = 0, Ai = 7r/2 only. Sometimes 
it happens that = for recorded data. In this case the 
data are phase insensitive yielding 1^=0 and the posterior 
phase distribution is homogeneous. 

Notice that this approach is well known in optics as phase 
estimator of discrete Fourier transformation (DFT) p^ . 
Define for a discrete signal Fj its DFT as 



X{m) 



1 J2 Fjexpi-i^njm/N). 



(20) 



As follows from the comparison of relations (pO|), ( p^ ) and 
( p^ , the visibility and phase of the generalized NFM treat- 
ment corresponds to the modulus and argument of the com- 
plex coefficient X{1) = R 



= arg{X(l)}, /^ = |X(1)|. 



(21) 



The discrete signal corresponds to the difference of regis- 
tered discrete scans of interference fringes created by chang- 
ing the auxiliary shift 



(22) 



The frequency m of the interference fringe in (21) is one 
since the fringe has been scanned only once (^4|). 

Thegeneralization NFM scheme given by phase estima- 
tion (|l^ is not optimal provided that the detected signal is 
Poissonian. In this case the likelihood function reads 



7V-1 
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(23) 



where the mean output intensities /" are given by ( |l3| ) and 
n'j represent the number of detected neutrons. Unfortu- 
nately, the explicit relations analogous to ([lO|-[T^) cannot 
be found analytically. The analysis must therefore be car- 
ried out numerically. However, various limiting and asymp- 
totic cases can be discussed analytically as shown in |18]. 
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Figure 1: Outline of the experimental setup. 



DFT estimation is the best estimation available provided a 
large number of particles is detected. No difference between 
Gaussian and Poissonian cases can be therefore expected 
in the regime of high intensities Io,Ih 3> 1. In the op- 
posite case of low output intensities Io,Ih ^ 1, the most 
frequent samples are "no detection at all" and "one neu- 
tron detected" in some beam at some value of the auxiliary 
shift. Direct substitution of this samples in both the Pois- 
sonian and the Gaussian likelihood functions tends to the 
prediction with undefined phase with = 0. 

Similarly to the previous case of strong field, where both 
methods are equally good, in the low field limit both meth- 
ods are equally bad. Similarly in the limit of low visibil- 
ity ^ min{ Jo, //i}, the results obtained with the use of 
the true Poissonian statistics are virtually identical to those 
yielded by the DFT |jl^. No difference between the wave 
and quantum approach can be observed here. The difference 
between the predictions for the quantum and wave phase es- 
timation in a realistic experiment should become significant 
for visibility close to unity and average number of detected 
particle of order one. 



4 Experiment 



Our experiments were performed at the neutron interfer- 
onietry setup at the 250 kW TRIGA reactor in Vienna. 
Thermal neutrons which are emitted from the moderator of 
the reactor behave like statistically independent particles. 
Therefore the correct description of the counting statistics 
of the input beam and both output beams is a Poissonian 
distribution p7|. Fig. 1 shows the experimental arrange- 
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position of the auxiliary shifter AG 



Figure 2: (a)Detected interference fringe 1° as mean of 690 single measurements (o) and 68.3% error bars for numerical 
values 1° = 2.21 neutrons, = 6.33 neutrons, = 1.03 neutrons and 6 = 4.83 rad. A typical single detection 1° denoted 
by full squares is shown as an example, (b) Experimentally obtained AE{A9) denoted by full squares are compared with 
theoretical prediction denoted by corresponding mean values (o) and error bars for 690 samples. 



ment. The input beam is split by the perfect crystal in- 
terferometer into two partially coherent beams. One of the 
beams passes a phase plate (grey shaded region) which in- 
troduces an unknown phase shift which has to be estimated 
from the experimental data. Then both beams are pass- 
ing an auxiliary phase shifter which modulates the output 
intensities at the detectors o and h. At the two output 
ports BF3 gas detectors enable single neutron counting with 
nearly 100 per cent efficiency. The very low intensities at 
the outgoing beams (1 neutron per second) allow a com- 
fortable electronical separation of the detector pulses. The 
mean number of collected neutrons is a linear function of 
the counting time which enables an adjustment of the de- 
sired intensities by proper selection of the counting time. 
In neutron interferometry an auxiliary phase shifter can be 
rotated in several discrete positions denoted by indices j in 
the intensity equation. Unique phase estimation is achieved 
even when other parameters of the setup (e.g. the mean 
intensities, the visibility and the frequency of the oscillation 
pattern) are unknown. In our experiment 8 equidistant po- 
sitions of the phase shifter were used for generation of the 
intensity modulation. 

To compare the efficiency of the NFM phase predictions 
with the optimum Poissonian MaxLik estimation the follow- 
ing procedure has been chosen. Each sample of data con- 
sisting of the number of neutrons counted in beams o, h in 
eight positions Aq, . . . , Ay (|lj), was processed using NFM 
formula (|l7| ) resulting in phase prediction Onfm- The rela- 
tive frequency fg{A6) characterizes how many times the es- 
timated phase Onfm faUs into the chosen phase window A9 
(confidence interval) around the true phase shift. The same 
procedure was repeated for phase prediction based on nu- 
merical maximization of the Poissonian likelihood function 
(H) §, I yielding the relative frequency of "hits" fp{Ae). 
The quantity 

AE{A0) = fp{Ae) - fg{Ae) (24) 



represents the difference in efficiency of the quantum and 
wave phase estimations for the given phase window A6. If 
this quantity is found to be positive, it means that the Max- 
Lik estimation is better than its Gaussian counterpart (sim- 
ply because more estimates of the phase shift fall into the 
chosen angular window if the former procedure is followed) . 
If, on the other hand, this quantity is not different from zero 
(in a statistically significant way) the two data evaluation 
procedures are statistically equivalent and no discrimination 
is possible. 

The result of the data analysis is shown in Figs. 2-4. Each 
figure consists of two different parts. The left panels show 
the detected (or simulated) data. The right panels then 
provide the interpretation of the corresponding data. Anal- 
ysis of the experimental data is summarized in Fig. 2. In 
addition to experiment, two Monte-Carlo simulations have 
been performed simulating experimental conditions under 
which the difference between Poissonian and Gaussian pre- 
dictions is negligible. The result of the simulations is shown 
in Fig. 3. Finally, the possibility to estimate several param- 
eters simultaneously is illustrated in Fig. 4. The difference 
AE was calculated using 690 experimental samples mea- 
sured in experiment with average beam intensities I°=2.2 
neutrons, /''=6.3 neutrons and visibility normalized with 
respect to o-beam being about 47%. As already explained 
in Section 2, this represents a critical situation, because even 
though there are less than 10 counts in each experimental 
run, one is nonetheless trying to get an estimation of the 
value of the phase shift. The experimental values of AE 
are depicted in Fig. 2b by full squares. For comparison, a 
theoretical prediction corresponding to the same values of 
parameters as in the real experiment were simulated in the 
Monte-Carlo experiment using 40000 samples. Open circles 
in the Fig. 2b show the corresponding mean values of the 
difference. However, since the experimental data are lim- 
ited due to the experimental conditions and available time 
to the relatively small number of 690 samples, the real data 
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Figure 3: (a) Interference fringe 1° (o) and error bars corresponding to computer simulation of the experiment with 
numerical values 1° = 2.21 neutrons, = 6.33 neutrons, = 0.258 neutrons and 9 = 4.83 rad. A typical sample 
denoted by full squares is shown as an example, (b) A£'(A^) corresponding to the same computer simulation is evaluated. 
Mean values (o) and error bars corresponding to 690 samples are given. The interpretation of panels (c) and (d) is the 
same as for (a) and (b), respectively. Numerical values are 1° = 0.551 neutrons, = 1.582 neutrons, = 0.258 neutrons 
and e = 4.83 rad. 



are fluctuating around the mean values. Statistical sig- 
nificance of the experimental results is demonstrated again 
using Monte-Carlo simulations. Other 20 simulations has 
been done, each of them with 690 samples The variance of 
the ensemble {Aii'j}, is shown in the Fig. 2b as "error bars" 
for each phase window. 

A significant difference between the effectiveness of clas- 
sical and optimum treatments is apparent in Fig. 2. The 
optimum estimation provides an improvement in fitting of 
the phase shift and the difference is beyond the estimation 
error, approximately 2.5 standard deviations in the opti- 
mum case. Obviously, no better performance of the MaxLik 
method can be expected for large values of the phase win- 
dow (any sensible statistical method would yield rather 
reasonable results). Likewise, no real improvement over the 
Gaussian estimate can be expected when A^ is close to zero, 
because too few data would then be accepted. The existence 
of a "best choice" for the phase window is therefore in itself 
an interesting feature of the method we propose. However 
notice that generalized NFM scheme (or equivalently DFT 
phase estimator) fits the phase shift quite well. The most 
pronounced difference is about 6% in the window of width 
about 1.256 rad. For example it means that the Gaussian 
fitting procedure hits this window 442 times, whereas the 
Poissonian one 484 times from a total number of 690 events. 
The experimental difference in the score is therefore 42 i.e. 
6% in favour of the latter method. This difference is a ran- 
dom number and theory predicts its value as 41 ± 12. The 



observed difference is therefore not large, yet statistically 
significant. 

In most cases, however, the difference between discrete 
and continuous nature of the signal is subtle enough to be 
hidden in various imperfections of the experimental setup. 
For illustration, a computer simulation of two experimental 
setups has been performed similarly to the above-mentioned 
case. A data analysis of the Monte-Carlo experiment with 
visibility as small as 1/4 of the visibility in Fig. 2 is shown 
in Fig. 3b. No statistically significant discrimination be- 
tween the Gaussian and Poissonian methods is possible in 
this case. The result of a simulation of an experiment with 
4-times less energy of interfering beams is shown in Fig. 3d. 
Also in this case, no discrimination is possible. In spite of 
the apparent interference patterns the classical description 
of phase can be fully justified both in the cases of low vis- 
ibility and low intensity. Particularly the latter case may 
appear as counterintuitive, since measurement with small 
number of particles is traditionally considered as a domain 
of quantum physics! 

Unlike the case of phase, the MaxLik estimation of the 
visibility is strongly biased in the case of small intensi- 
ties. There is simple explanation for this behaviour. In 
low-intensity regime the character of individual detected 
samples is determined rather by fluctuations then by ac- 
tual parameters of the experimental setup as seen in Fig. 3c 
for example. The MaxLik estimation of the visibility fits 
this fiuctuations and, as a consequence, the estimation is 
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Figure 4; (a) Interference fringe 1° (o) and error bars corresponding to a computer simulation of the experiment with 
numerical values 1° = 22.1 neutrons, = 63.3 neutrons, = 10.3 neutrons and 9 = 4.83 rad. A typical sample is 
denoted by full squares as an example, (b) Histogram of estimated visibilities normalized with respect to obeam obtained 
in the same simulation. The true normalized visibility is 46.7%. 



biased. This is particularly obvious in the case 
when interference pattern disappears, but MaxLik fitting 
yields a visibility of about 40% due to fluctuations. Never- 
theless, except the case of very low intensities, the estimated 
visibility is meaningful as demonstrated in Fig. 4, where a 
histogram of the estimated visibilities normalized with re- 
spect to o-beam is shown as result of a simulation with 10 
times stronger output beams compared to the output beams 
used in the experiment. It is apparent that the estimated 
visibilities are distributed with Gaussian-like shape around 
the value of the true visibility. 

5 Conclusion 

A statistically motivated analysis of neutron interferome- 
try provides a correction to the previously introduced op- 
erational quantum phase concept. Since the standard ap- 
proach has been universally derived, without considering 
the statistics of the interfering fields, it cannot be optimal. 
This additional knowledge may be used for improved pre- 
dictions and testing. This scheme provides therefore a sta- 
tistically motivated evaluation of the whole interferometric 
system. Instead of the question of the wave theory: "How 
precisely can the interference fringes be distinguished?" a 
more sophisticated question is here formulated as: "What 
statistical properties can be recognized from an interference 
pattern?" In particular, the experiment performed with neu- 
trons demonstrated a measurable improvement of phase fit- 
ting for discrete Poissonian signals. 
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